Non-equilibrium effects in a Josephson junction coupled to a precessing spin 
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We present a theoretical study of a Josephson junction consisting of two s-wave superconducting 
leads coupled over a classical spin. When an external magnetic field is applied, the classical spin will 
precess with the Larmor frequency. This magnetically active interface results in a time-dependent 
boundary condition with different tunneling amplitudes for spin-up and spin-down quasiparticles and 
where the precession produces spin- flip scattering processes. We show that as a result, the Andreev 
states develop sidebands and a non-equilibrium population which depend on the precession frequency 
and the angle between the classical spin and the external magnetic field. The Andreev states lead to 
a steady-state Josephson current whose current-phase relation could be used for characterizing the 
precessing spin. In addition to the charge transport, a magnetization current is also generated. This 
spin current is time-dependent and its polarization axis rotates with the same precession frequency 
as the classical spin. 

PACS numbers: 



I. INTRODUCTION 

Recently, superconducting- ferromagnetic (SF) hybrid 
devices have received increased attention due to their 
potential as spintronics devices. In spintronics, the spin 
degree of freedom is employed to create new phenomena 
which could be used to create entirely new devices or 
be used in combination with conventional charge-based 
electronical'^. Information, e.g., can be stored in the 
magnetization direction of a small ferromagnet and its 
state can be read out by measuring a current through 
a nano-scaled contact determined by the magnetization 
direction. Nanomagnets such as single molecular mag- 
nets or magnetic nanoparticles may be suitable building 
blocks for such information storage^ ^. 

The interest in single molecular magnets been sparked 
by their appealingly long relaxation times at low 
temperatures^ and experimental breakthroughs in con- 
tacting molecules to both superconducting and normal 
leads has made molecular spintronics a growing field of 
research. Transport measurements of molecular mag- 
nets in normal junctions have been made as a means 
to characterize the magnetic states^ ^. Contacting o f Ceo 
molecules^, metallofullerenes^ and carbon nanotube d"*^^ * * 
to superconducting leads has also been demonstrated. 
In addition, single molecular magnets ha ve bee n sug- 
gested for quantum computing application d^^l^^l due to 
their long relaxation times. 

Currents are not only used to read out the state of a 
magnet but are also used to control the magnetization 
direction. Spin-polarized currents carry angular momen- 
tum. However, a spin current is not a conserved quantity 
in a ferromagnet and a spin current oriented in such a 
way that its direction is perpendicular to the interface 
plane between the ferromagnetic layer and the leads may 
lose some of its spin-angular momentum. The angular 



momentum lost by spin-polarized electrons transported 
through a ferromagnet is transferred to the ferromagnet. 
This transfer of angular momentum generates a torque 
acting on the feromagnet's magnetization direction. This 
spin-transfer torque mediated by electrical currents was 
theoretically investigated by Slonczewski^ and Bergei'^ 
who worked out a description for ferromagnet-normal 
metal (FN) multilayer structures and showed that spin- 
transfer torques can lead to precession as well as reversal 
of the magnetization direction. These theoretical predic- 
tions were experimentally verified by Tsoi^^ and Myer4^. 
Nonequilbrium magnetization dynamics^^ and spin- 
transfer torques^^ in FNF trilayers coupled to supercon- 
ducting leads have also been studied. 

In this paper, we study the coupling between the mag- 
netization dynamics of a nanomagnet or single molecular 
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FIG. 1: (Color online) Two superconducting leads are coupled 
over a magnetic spin, S{t). The spin precesses with the Lar- 
mor frequency ul due to an applied external magnetic field 
H applied at an angle i9. Quasiparticles tunneling between 
superconductors interact with the precessing spin via an ex- 
change coupling resulting in a dynamical inverse proximity 
effect existing within the superconducting coherence length 
^0 of the junction interface. 
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magnet and Josephson currents through a nano-scaled 
junction. We consider two superconducting leads cou- 
pled over a nanomagnet, consisting of e.g. a molecular 
magnet or a magnetic nanoparticle as shown in figure [T] 
The spins of the magnetic molecule or nanoparticle are 
assumed to be held parallel to each other resulting in 
a uniform magnetization which can be represented by a 
macrospin^^. If an external magnetic field is applied, the 
spin of the nanomagnet starts to precess with the Lar- 
mor frequency. This dynamics changes howeve r whe n it 
is coupled to conduction electrons in the lead d^^ * ^ ^ . As 
our starting point, we take the model by Zhu et al^^^ 
and extend it to include arbitrary tunneling strengths 
leading to a modified quasiparticle spectrum displaying 
Andreev levels for energies within the superconducting 
gap, l^l < A(T)^^^^. In reference [26^, we focused on 
the dc Josephson charge current, while here, we focus on 
the coupling between the dynamics of the Andreev levels 
and the transport properties. The coupling of the two 
superconducting leads over the precessing spin produces 
an ac spin Josephson current. The difference between 
the spin currents on the left and right sides of the in- 
terfaces produces a spin-transfer torque, r = j£ — Jr^ 
shifting the precession frequency of the rotating spin. At 
finite temperatures, there is also a spin current carried by 
quasiparticles generating a damping of the magnetization 
dynamics, the so-called Gilbert damping^^ A transi- 
tion of the leads from the normal state into the super- 
conducting state reduces the Gilbert damping^ since the 
number of quasiparticles is suppressed for temperatures 
T < tJ^. This interplay between the Josephson effect 
and a single spin may be used for read-out of quantum- 
spin states^^ or for manipulation of Andreev levels in the 
junction^. 

Furthermore, we find that the ac Josephson spin cur- 
rent is a result of superconducting spin-triplet correla- 
tions induced by the spin precession. The appearance 
of superconducting spin-triplet correlations in SFS junc- 
tions has been used to explain the observation of a long- 
range proximity effect in a number of experiments^^^^. 
Keizer et al. observed a supercurrent through a junc- 
tion consisting of conventional s-wave superconductors 
coupled over a layer of the half-metallic ferromagnet 
Cr02 much thicker than the decay length of the super- 
conducting spin-singlet correlation^^. Various mecha- 
nisms for converting the spin-singlet correlations of the 
superconducting leads into spin-triplet correlations that 
may survive within a ferromagnetic layer have been 
suggestecP^^. Bergeret et al. showed that a local inho- 
mogeneous magnetization direction at the SF interface is 
sufficient to generate spin-triplet conversiori^. In refer- 
ence [39], it was suggested that the spin-singlet to spin- 
triplet conversion is due to interface regions with mis- 
aligned averaged magnetic moments breaking the spin- 
rotation symmetry of the junction producing spin mix- 
ing as well as spin-flip processes. A similar trilayer 
structure with noncollinear magnetizations resulting in 
a long-range triplet proximity effect was proposed by 



Houzet and Buzdin^. Taking into account the impor- 
tance of interface composition, Khaire et alW^ devised 
SFS junctions consisting of conventional superconductors 
and Cr02 in which they had inserted weakly ferromag- 
netic layers between the superconductors and the half 
metal to produce interface layers with misaligned mag- 
netization directions. A long-range proximity effect was 
observed in junctions containing the interface layers, but 
not in junctions without. Conflrmation of Keizer's re- 
sults were made by Wang et al.^^ who measured a su- 
percurrent through a crystalline Co nanowire. The Co 
nanowire was a single crystal, but the contacting proce- 
dure was likely to cause defects at the SF interfaces and 
the inhomogeneous magnetic moments needed to create 
the spin-triplet correlations. Other experimental verifl- 
cations of long-range proximity effects includes Holmium 
(Ho) wires contacted to conventional superconductor^^. 
Ho has a conical ferromagnetic structure whose magne- 
tization rotates like a helix along the c axis. The ap- 
pearance of spin-triplet correlations in such junctions and 
their effect on the long-range proximity effect^^ and 
spin current have also been studied theoretically. In 
the present problem, the magnetization direction varies 
in time rather than in space giving rise to time-dependent 
Andreev level dynamics and a dynamical inverse prox- 
imity effect in the form of induced time-dependent spin- 
triplet correlations. Houzet^^ studied a related problem 
in which a Josephson junction consisting a ferromagnetic 
layer with a precessing magnetization placed between two 
diffusive superconductors was predicted to display a long- 
range triplet proximity effect. 

We formulate the problem of two superconduct- 
ing leads coupled over a nanomagnet in terms of 
nonequilibrium Green's functions. The quasiclassi- 
cal theory of superconductivity is based on Landau's 
Fermi liquid theor^H^^ and is applicable to both 
superconducting^^^t^ and superfluicP^ phenomena as 
well as inhomogeneous superconductors and nonequilib- 
rium situations. The quasiclassical theory gives a macro- 
scopic description where microscopic details are entered 
as phenomenological parameters'^. Basically, it is an ex- 
pansion in a small parameter kBT/Ep^ where is the 
Fermi energy, and is suitable for weakly perturbed super- 
conductors. The perturbations should be weak compared 
to the Fermi energy, V <C Ep^ and of low frequency, 
fvjo <C Ep- Interfaces and surfaces in superconducting 
heterostructures or point contacts, on the other hand, are 
strong, localized perturbations with strengths compara- 
ble to the Fermi surface energ}''^. Within quasiclassical 
theory, interfaces are handled by formulation of boundary 
conditions which usually have been expressed as scatter- 
ing problems, being able to treat spin-independent^^^^^^ 
as well as spin-dependent, or spin-active, interfaces'^''^. 
In many problems, in particular when an explicit time 
dependenc e appea rs, the T-matrix formulation is more 
convenient ^^'^^'^. This formulation is also well suited 
for studying interfaces with different numbers of trajec- 
tories on either side as is the case for normal metal/half 
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metal interfaces^^^^. The two methods have proved to be 
equivalent and may be applied both in the limit of clean 
and in the limit of diffuse superconductor^^. In the lat- 
ter case, the boundary conditions coincides with those 
of Kuprianov and Lukiche\^^ and of Nazaro\^^. For a 
recent review of quasiclassical theory we refer the reader 
to reference [67 . In the present problem, the dynamics 
of the nanomagnet constitutes a time-dependent spin- 
active boundary condition for the two superconductors 
which we solve using the T- matrix formulation. First, 
the transport equations are solved separately to find the 
classical trajectories for each lead. Then the T-matrix 
describing the scattering between the leads is used to 
connect the trajectories across the time-dependent spin- 
active interface. 

We start by outlining the T-matrix formulation appli- 
cable to scattering via the precessing magnetic moment 
in Sec. II. We show that the boundary condition can be 
solved both in the laboratory frame and a rotating frame. 
In the latter solution, the explicit time dependence is re- 
moved by a transformation to a rotating frame rendering 
this approach suitable for efficient numerical implemen- 
tations for studies of transport properties. However, the 
solution comes at the cost of introducing an energy shift 
of the chemical potentials for the spin-up and spin-down 
bands. The laboratory frame approach is, on the other 
hand, suitable for studying modifications to the super- 
conducting state although the explicit time dependence 
increases the complexity of the solution. In Sec. Ill A, 
we review the results for the Josephson charge current in 
reference [26] in terms of the laboratory frame descrip- 
tion. The spin currents are described in Sec. Ill B, which 
is followed in Sec. Ill C by the induced time-dependent 
spin-triplet correlations and Andreev-level dynamics giv- 
ing rise to the spin currents. In Sec. Ill D, we discuss 
the back-action of the scattering processes on the mag- 
netization dynamics while the magnetization induced in 
the leads is discussed in section III E. In Sec. IV, we 
conclude with a summary of our results. 



II. MODEL 

We consider two superconductors forming a Josephson 
junction over a nanomagnet. The nanomagnet may ei- 
ther be magnetic nanoparticle or a single-molecule mag- 
net and we will assume that contact between the leads 
and the nanomagnet is made up of a few single quantum 
channels. The magnetization of the nanomagnet is put 
in precession and the resulting contact will constitute a 
time-dependent spin- active interface (see figure [T]). The 
nanomagnet together with the two superconducting leads 
are described by the total Hamiltoniarl^^^ 

n = nieads^nB^nT- (i) 



The left (L) and right (R) leads are s-wave superconduc- 
tors described by the BCS Hamiltonian 

a = L,R a = L,R 

(2) 

where the dispersion, (f/. = Ti?k'^ /2m — /i, and the chem- 
ical potential, /i, are assumed to be the same for both 
leads. The order parameter of the leads is assumed tem- 
perature dependent, Aq, = A(T)e^*^/^. Here Lp is the 
relative superconducting phase difference over the junc- 
tion which we treat as a static variable that is tunable. 
The nanomagnet is subjected to an external magnetic 
field modeled as an effective field H acting on the nano- 
magnet's magnetic moment, ijl. Included in this effective 
field are also any r.f. fields to maintain precession, crystal 
anisotropy fields and demagnetization effects. The mag- 
netic moment of the nanomagnet is viewed as a single 
spin, or macrospin, which we will treat as a classical en- 
tity. This macrospin is related to the magnetic moment 
hy II = —jS where 7 is the gyromagnetic ratio. The 
spin and the effective magnetic field couple via a Zeeman 
term, 

= -jS • H. (3) 

If the effective field is applied at an angle, i?, relative to 
the spin, a torque is produced that brings the classical 
spin into precession around the direction of the effective 
field. This precession generated by the tilt angle occurs 
with the Larmor frequency ujl = 7^, where the mag- 
nitude of the external field is H = \H\. Here, we take 
the direction of the effective field to be along the axis 
and the angle 1^ is defined diS H ■ S{t) = HS cos 1^. The 
Larmor precession is captured by the equation of motion 

f = -,SxH (4) 

where the right-hand side is the torque produced by the 
effective field. The magnitude of the spin is constant, 
S = \S{t)\ = const. Consequently, the path traced out 
by the spin, S{t) = Ses{t), has a time dependence which 
lies totally in the direction of the spin, es{t). In the 
absence of any other torques than the effective field, the 
direction of the precessing spin may be written as 

^sit) = {cos{ujLt) sini^e^^ + sin(co'L^) siwd Cy + cosi^e^). 

(5) 

Tunneling quasiparticles have the possibility to tunnel 
directly between the leads with hopping amplitude Vq or 
interact with the precessing spin via an exchange cou- 
pling of strength Vs. These processes are described by 
the tunneling Hamiltonian 

k(j\k' a' 

(6) 

where Vka;k'^' = {VoS^^, + Vs{S{t) ■ (T)^„,)S{k - k') and 
cr = {crx,cry,crz) with cr^ being the Pauli matrices. The 
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spin-dependent hopping amplitude is 

VsS{t)-(T = 1/55' ( COS ^ cr^ + sin ^e-^^^^^^cr^) (7) 

where the first term is a spin-conserving part with dif- 
ferent hopping amplitudes for spin-up and spin-down 
quasiparticles. The second term is time-dependent and 
describes processes where quasiparticles flip their spins 
while exchanging energy ujl with the rotating spin. The 
junction's transport properties as well as the modifica- 
tions to the superconducting states of the leads depend 
on the hopping amplitudes Vq^Vs as well as the the su- 
perconducting phase difference, (p, the precession of the 
spin, S{t), and the tilt angle, 1^^ between the effective 
field and the precessing spin. 



A. Approach 

We formulate the problem using nonequilibrium 
Green's functions in the quasiclassical approximation fol- 
lowing references [58 64 . The tunneling Hamiltonian (|6| 
provides a time-dependent and spin-active boundary con- 
dition for the quasiclassical Green's function and is solved 
by a T-matrix equation^^**^^ . A quasiclassical Green's 
function is a propagator describing quasiparticles moving 
along classical trajectories defined by the Fermi velocity 
vf = vf{Pf) a given quasiparticle momentum on the 
Fermi surface, p^. The information of a quasiclassical 
Green's function is contained in the the object 



g{PF,R;e,t) 



g^{pF,R',e,t) 



(8) 



where R is the spatial coordinate, e is the quasiparti- 
cle energy relative to the chemical potential /i and t is 
time. The g{pp^ R;e^t) propagator is a 8 x 8 matrix in 
the combined Keldysh-Nambu-spin space. The "check" 
denotes a 2 x 2 matrix in Keldysh space where the com- 
ponents are retarded (i?), advanced {A) and Keldysh (K) 
Green's functions while the "hat" indicates a 2 x 2 ma- 
trix in Nambu, or particle- hole, space which is further 
parameterized using the Pauli spin matrices {(Jx^o-y^az) 
(see reference [50] for details). The matrix components in 
the combined Nambu-spin space are conveniently divided 
into spin scalar (s) and spin vector (t) parts. 



.,x 



{X = R^A^K) where and are the spin-singlet 
and spin-triplet components of the anomalous Green's 
functions (similarly for g^ and g^). The Green's 
function g(pj^, R;e,t) obeys a Boltzman-like transport 



equation 



4714491 



ihvF-Vg + [ehl - U,g]o = 0, 



where are Pauli matrices in Nambu space and i-L in- 
cludes self-energies like the superconducting order pa- 
rameter, A = Al, and impurity contributions, T^imp^ as 



well as any external fields, hext- The "o" product repre- 
sents a convolution over common time arguments com- 
bined with a matrix muliplicatiorP^'^. The transport 
equation is complemented by the Eilenberger normaliza- 
tion condition, 



gog = -TT^i, 



(11) 



and a set of self-consistency equations such as the one for 
the superconducting order parameter. 



A(i^,t) = A 



de 



where is an average over the Fermi surface, A is the 

pairing-interaction strength and Sc is the cut-off energy 
which may be eliminated by making use of the critical 
temperature, Tc. 



The transport equation (10) is solved separately for 



each lead, treating the interface as an impenetrable sur- 
face where quasiparticles are perfectly reflected. This 
hard wall boundary condition leads to a solution, 
for each semi-infinite lead, a = L^R. The propagators 
are then connected across the interface by the tunneling 
Hamiltonian, Ht, whose effects can be incorporated via 
a quasiclassical t-matrix equation as 

ta{t,t') = t^{t,t') + [t^ogloQ{t,t'). (13) 

The hopping elements of the tunneling Hamiltonian enter 
the t-matrix equation via a matrix TL{t^t') defined as 



rL{t,t') = [vog%ov]{t,t') 



(14) 



for the left side of the interface and the right-side matrix 
is obtained by interchanging L and R. The time 
dependence in the current problem enters through the 
hopping element v{t)^ which in particle- hole space has 
the form 



Vo 



vses{t)-(T 
Vo 





vs':ry{e.s{t)-cr)ay 



(15) 



and v{t) = v{t)l in Keldysh space. Here, the hop- 
ping elements Vq^Vs in equation ([6| have been replaced 
with their Fermi-surface limits, vq = ttNfVo and vs = 
ttNfSVs with Nf being the normal density of states at 
Ef' Note that for the junction studied, the hopping ele- 
ments have the symmetries VLR{t) = ^lr{'^) = ^Rhi't) = 
v(t). The t matrices (13) are used to calculate the full 
quasiclassical propagators which depending on if their 
trajectories lead up to or away from the interface are 
divided into "incoming" {g^) and "outgoing" {g^) propa- 
gators, given by 



(10) 5r(t,i')=5a(i,i') + [(ff°±»7ri)ot„o(50Ti7ri)](t,t') (16) 



where a = L^R and the upper and lower signs, ± and ^, 
refer to the incoming and outgoing propagators, respec- 
tively. 
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The self-energy fields, such as the order parameter, de- 
pend on the full propagators and should in principle be 
calculated self-consistenly taking into account the inter- 
face scattering. However, we assume in this study that 
the area of the point contact, is small compared to tt^q 
where = ^|t^F|/27rTc is the superconducting coherence 
length. As a result, the superconducting state does not 
change considerably and the order parameter. A, and 
other possible self-energies in the leads do not have to be 
recalculated^^. We will also assume that the supercon- 
ducting phase changes abruptly over the contact. 

The use of equations ( 1Q|11 ) together with the bound- 
ary condition (16) allows for calculation of the transport 
properties the junction. The charge and spin currents 
are given by an average over the Fermi-surface momen- 
tum directions of the full propagators. Here, this average 
amounts to a difference between incoming and outgoing 
propagators and the charge, j^, and spin, j^, currents 
evaluated in lead a for a single conduction channel are 

m = ^1V[f3(5L'<(e,t) -5°'<(£,t))] (17) 

and 

Jait) = \j |^Tr[f3<T(5^<(£,t) -5°'<(e,t))], (18) 

where & = diag(t7, — a^ytrcr^). The Green's functions in 
the above expressions are the lesser propagators defined 



as g 



domain to energy space where it becomes an algebraic 
equation. 



The propagators g^{£^£') = g^{£)5{£ — £') have the fol- 
lowing Nambu-spin structure. 



TT 



.R 



A(T, ±v)iay 



QR \iayA*{T,±ip) -e' 



where fi^ = y'|A(T)|2 - (£-R)2, e« = e + io+, 

=T3[5°'^(£)]^f3, (20) 

and 

^S'^(^) = (^) - (^)) tanh(e/2T). 

The gap, A(T,±(/p) = A(T)e±^'^/^ is both tempera- 
ture and phase dependent and the "+"("—") sign of 
the phase dependence refers to lead R{L). Furthermore, 
tanh(£:/2T) is the quasiparticle occupation function set- 
ting the two superconducting leads in thermal equilib- 
rium with each other. The t matrices in energy space 
are a sum of t matrices whose energies differ hy ujl and 
satisfy the relation 



ia{£,£') = ^ta{£,£ ^nUJL)S{£ - 



' + ncjL), (21) 



B. Solving the time-dependent boundary condition which is equivalent to a time dependence of 



The time-dependent boundary conditions can be 
solved in two different ways. In the first procedure, 
the boundary conditions are solved in the laboratory 
frame in which the time dependence is preserved and 
is manifested as frequency shifts in a difference equa- 
tion. The treatment is similar to that of dc-voltage bi- 
ased SIS junctions^o ^8 59 63 711MI. The second approach 
involves removing the explicit time dependence of v{t) by 
a transformation to a rotating frame (see reference [26 
for details). This procedure is numerically more efficient 
but the transformation, however, introduces an exchange 
field shifting the chemical potentials of the spin-up and 
spin-down bands in the leads making the first approach 
more suitable for studying changes in the superconduct- 
ing state in vicinity of the junction due to quasiparticle 
tunneling via the precessing spin. Below we describe both 
within a quasiclassical framework. 



1. Laboratory frame 



To solve the boundary condition ( 16 ) dependent on the 



matrix (14) in the laboratory frame, it is more convenient 



-inojLt 



I 



d£ 
2^' 



-i£{t-t')l 



(22) 

The Ansatz above and the assumption that the leads are 
in equilibrium so that their respective Green's functions 
may be written as g^(t^t') = g^{t ~ t') lets one evaluate 
the coefficient matrices ^^^(e,^') and ra{£^£ )ga{^ ) 



equation (19) and subsequently solve t he result ing differ- 
ence equation in terms of ia{£ + ncj^ , ^)^^^'. This is a 
quite general procedure capable of handling diverse forms 
of self-energy fields H. However, the matrix coefficients 
in equation (19) must be evaluated for each particular 
kind of junction and lead state. The properties of these 
matrices then determine the specific t-matrix difference 
equation and the solution strategy. 

For the present calculation certain simplifications can 
be made due to the spin independence of g^{£) given 
by equation (20) and the form of the hopping element 
in equation (15): the Keldysh-Nambu-spin matrices can 



to Fourier transform the t-matrix equation from the time 



be factorized in spin space into generalized diagonal ma- 
trices X^, spin-raising matrices X^, and spin- lowering 
matrices X^. These matrices are still Keldysh-Nambu 
matrices but they have the algebraic properties of spin 
matrices such as X^oY^ X^oY^ 0, X^'^oF^'^ oc Z^, 



and o y^'^ oc Z^'^. A matrix factorized in this form 
may be shown to have the time dependence 

de 



-ie(t-t') 



27r 



(23) 



-iojLt 



XHe, 



Using the spin algebra, the t-matrix equation in energy 



space (equation ([19|) may be written as 

Af 



1 



^0 





1 



- 

^0 



-Bt 



1 






/ft 











(24) 



The coefficient matrices in equation (24) are functions 
of energy e and precession frequency ul and are straight 
forward to evaluate. Contrary to e.g. the case of a finite 
dc voltage, where the t-matrix equation is a difference 
equation solved by recursive methods, we have a matrix 
equation for t which can be solved by simple (numerical) 
inversion. Factorizing the propagators in equation (Il6| 
according to their spin structure results in 

ai^'^ie) =gl{e) + Mi±{^)i'i{e,ujL)Mt,^{e) (25) 
5l-^.°(e,t)=e---*M^,±(£+c^^)il(£,c.i)M„^^^(e)(26) 
5i.^°(e,t)=e+--*M„^,±(e-Wi)tt(e,Wi)M^,^(£)(27) 

with Mi±{e) = gl{e)±ml. 



2. Rotating frame 

The unitary transformation matrix for removing the 
time dependence of v{t) is 











(28) 



resulting in a transformation to a rotating frame of ref- 
erence with respect to the precessing spin S{t) in which 
the hopping element is 



- vses-cF 
^0 







v=U\t)v{t)U{t) = 

(29) 

The direction of the precessing spin is now static 
in this rotating frame, = cosi^e^ + sin'^^e^;, but the 
hopping element v is still spin active with different hop- 
ping amplitudes for spin-up and spin-down quasiparti- 
cles, ^0 i ^s-cosi?. The hopping element also contains 
a spin- flip term vs sin ^ scattering between the two spin 
bands. 

Next, we apply the unitary transformation to the prop- 
agator go,{t^t') and obtain 



(30) 



de de' 



g^{e,e')U{t'). 
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FIG. 2: (Color online) (a) Schematic illustration of Andreev 
scattering in a Josephson junction. Electron- like quasipar- 
ticles (blue) impinging on the right superconductor are An- 
dreev retroreflected (thick arrows) into hole-like quasiparticles 
(white) at the superconductor surface. The hole-like quasi- 
particles are transmitted (thin arrows) to the opposite (left) 
superconductor, where they are Andreev retroreflected into 
electron-like quasiparticles. This Andreev scattering process 
is phase coherent and results in the formation of quasipar- 
ticle states within the superconducting gap, Andreev levels, 
when the quasiparticles interfere constructively along close 
loops, (b) Scattering processes in a Josephson junction cou- 
pled coupled over a classical spin precessing with frequency 
(jJl- a tunneling quasiparticle has the possibility to tunnel 
across the junction with preserved spin and energy, or having 
its spin flipped by interaction with the classical spin while 
gaining or losing energy ul- Andreev reflection occurs at the 
junction interfaces. 



The transformation of the propagators into the rotating 
frame introduces spin-dependent energy shifts, e ^ e ^ 
displayed in the transformed propagators, e.g. as 



me) 



(31) 



gSie- 



2 / 
2 / 



gSie- 



2 / 
2 / 











for the retarded Green's function. The advanced and 
Keldysh propagators are similarly transformed. 

Applying U{t) to the transport equation, equation 



(10), similarly leads to a shift of the energies of the spin- 



up and spin-down quasiparticles. The resulting transport 
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equation is 

ihVp-Vga + [{STs + heff (7^)1 - Ua^QaV = 0, (32) 

where the energy shifts are captured by the appearance 
of an effective magnetic field in the leads, h^ff = ^e^. 
Equation (32) is time independent as long as the self- 
energy fields in do not contain terms which are off- 
diagonal in spin space corresponding to spin-flip scatter- 
ing or equal-spin paring order parameters. For ballistic, 
or clean, s-wave superconducting leads. Ha = ^(^)q; 
and equation (|32|) is time independent. The propaga- 



tors ga{t^f) obey e quati on (32) for a specular-scattering 
boundary conditiorP^^. Introduction of the tunneling 
processes lead to a boundary condition problem which is 
still time independent in this rotating frame and the " o" 
product in equation ( 16 ) reduces to a matrix multiplica- 
tion in energy space. 

The precessing spin introduces spin mixing of the two 
spin bands in such a way that spin- up (-down) quasi- 
particles are scattered or injected into the spin-down (- 
up) band. This non-equilibrium spin injection, however, 
leaves the charge current, equation (17), time indepen- 



dent. The charge current's steady- state solution follows 
directly from the spin-independent trace over particles 
and holes in equation (17) and the fact that the trans- 



formation matrix U{t) leaves the diagonal terms of the 
propagators ga^'^ time independent. The spin current, 
on the other hand, has components which include the 
off-diagonal elements of the propagators ga^'^ . The off- 
diagonal elements, which are proportional to the Pauli 
spin matrices ax and ay, are time dependent since these 
are affected by the time dependence of the transforma- 
tion matrix U{t). Thus, for a finite tilt angle, ^ ^ 0, the 
spin current may be time dependent. 



III. RESULTS 

The precessing spin introduces a new energy scale 
and the superconducting correlation functions can be ex- 
pected to be modified due to the scattering processes the 
precessing spin gives rise to. The starting point is the 
s-wave superconducting leads and their the spin-singlet 
pairing amplitudes, fs ((V^/ctV^-fc; - ^l^-kii^kt))' When 
electron- and hole-like quasiparticles interfere construc- 
tively, sharp states inside the superconducting gap called 
Andreev states are formed. In regular Josephson junc- 
tions without magnetically active interfaces, these An- 
dreev states come in degenerate pairs V^^^), which can 
be described by the spinors V^^^) = ('^/'^(i), V^j^^^^)"^. Each 
member of the spinor, V^^^), is subjected to transmission 
yLR{RL) and Andreev retroreflection ^i?(L) and the An- 
dreev bound states are formed when the processes lead 
to constructive interference along closed loops schemat- 
ically illustrated as ^l)L^^^){k,£) ^ V^i,^(;)(/c, e) ^ 



schematic picture of the scattering processes is shown in 
the upper panel of figure [2] The tunneling processes de- 
scribed above are captured by the hopping element 



Vo 



- vs COS i^az 






vq + Vs COS i9az 



(33) 



which lets quasiparticles tunnel across the junction with 
their spin directions and energies unaffected but with 
different hopping amplitude for spin- up and spin-down. 
This scattering behavior is captured by the matrix o 
gn o (see below). The spin- flip part of the hopping 
element, 



Mi) 



Vs sin 'i^cr+(_ 






Vs sin i^a_ 



•(+) 



(34) 



flips a spin-down (-up) quasiparticle into a spin-up (- 
down) quasiparticle while changing its energy by 
ujl{—^l)' Here, we have defined cr+(-) = ^{o'x i 
iay). This kind of tunneling creates spin- flip pro- 



£_(£+)) ^ 

-k-e) ^ 



cesses, e.g. ipL^(i.)(k,e 

i^uni)^-k-{-k+),-e-{-e+)) "-^ 
'0Lt(i)(^5 which is a process captured by the elements 

Focusing on the left side of the interface and pa- 



rameterizing the matrix F^ according to equation (23), 
leads one to conclude that the spin-preserving compo- 
nent of Fi, is a combination of the two processes de- 



V^og^ov^ 



V^og^ov^ . 



scribed above, Fj = v'^og'^ov 
There are also mixed tunneling processes where tunneling 
with and without spin flip are combined into the terms 



V^og^cyv^ 



-v'^og^ov^ and f;^ - 



v^og^av^ - 



-V^og^av^. 



These two terms generate a net spin-flip for quasiparti- 
cles tunneling across the nanomagnet. 

The angle "3 between the spin S and the magnetic field 
H determines the amount of spin-flip scattering ranging 
from zero for parallel alignment to maximum in the case 
of ± H. The frequency ujl sets the amount of energy 
exchanged between a quasiparticle and the rotating spin 
during a spin- flip event as indicated in figure [2] (b). First, 
we will look at the consequences for the density of states 
and the charge current due to the scattering caused by 
the precessing spin. After that, we will take a closer look 
at the effects on the superconducting pair correlations 
before we turn to the spin scattering states and the spin 
current as well as their implications for the leads. 



A. Charge currents 

This section reviews some of the results presented in 
reference [26]. Here, however, the charge current results 
are described in the t-matrix formulation and are in- 
cluded for completeness. 
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function aPEIl 



■c 2e dej , . £j ^ 



(37) 



which, given the energy dispersion in equation (36), is 
evaluated to 



3a 



{Vo-Vs)Asmip 



1 - Do sin^ {ip/2) - Vs cos^ {ip/2) 



tanh ( ^ ) . 

(38) 



I showing that he critical current is reduced due to 
°^ the spin-flip scattering generated by the embedded 
nanomagnet^^. 

If only spin-independent tunneling is present and the 
■05 spin-dependent hopping strength vs 0, the spin- 
transmission coefficient Vs while Vq ^ V = 



FIG. 3: (Color online) Current-phase relation (a) and An- 
dreev levels (b) for a junction whose hopping elements are 
parameterized as vq = cos a and vs = sin a and a is var- 
ied from (solid black line) to 7r/2 (solid red line) in steps 
of 0.l7r/2. The charge current is plotted in units of eA/h 
for temperature T ^ 0. The charge-current kernel j^'^{s, cp) 
shows the population of the Andreev levels as well as the di- 
rection in which the populated Andreev levels carry current. 
In (c), j^'^ is plotted for a = 0. The lower Andreev level, 
which has energies below the Fermi surface, is populated and 
carries current in the positive direction (red) for phase dif- 
ferences < (p < TV and in the negative direction (blue) for 
TT < cp < 27T. In (d), j^'^ is similarly plotted for a = 7r/2. 
In this case, the populated Andreev level carries a negative 
current for < cp < tt and a positive one for tt < cp < 27t. 
The lower gap edge is indicated by the dashed black line. 



1. Static spin 

For a static spin with precession frequency ujl 
t-matrix equation ( 24 ) reduces to 



[1 - Ag] = F' 



0, the 



(35) 



which can easily be solved analytically since the spin- 
up and spin-down bands separate into two sets of equa- 
tions. Straight-forward calculations of the density of 
states show that Andreev levels form within the super- 
conducting gap and are located at energies given as 



£j = ±AW 1 — Vq sin 



2 ^ 



Vs cos^ 



(36) 



where we have defined Vq = 4'Uq/(1 + 2{vq + Vg) + ('^q 
v's)'^) and a spin-transmission coefficient Vs = 4'u|/(l 



The charge current density 



is related to the energy dispersion and the occupation 



,2\2 



which is the usual transparency of a 



Josephson junction. The Andreev levels are now sj 
±Aa/1 — Psin^ ^ and carry a charge current 



Ja 



VA sin (f 



1 -Psin^ {(p/2) 



tanh (^). 



(39) 



Both of these relations are shown as solid black lines in 
figures [3] (a) and (b). As can be seen in these figures, 
the junction is in a state corresponding to the junc- 
tion's energy being minimized for the phase difference 
(p = 0. Increasing the hopping strength of the spin- 
dependent tunneling until the spin-dependent tunneling 
dominates, vs > vq^ causes the junction to shift from be- 
ing in the state to being a tt state, as can be seen in fig- 
urejsja). When spin-dependent tunneling dominates, the 
junction's ground state is such that the coupled super- 
conductors have an internal phase shift of tt, as predicted 
in reference ^77j. Such tt states can also be observed in 
junctions where the spin- active barrier has been extended 
to a ferromagnetic region^^. In such junctions, the width 
of the ferromagnetic layer as well as the strength of the 
exchange field determine the transport properties. 

If, on the other hand, the spin-independent tunneling 
is decreased to 0, leaving only spin-dependent tunneling, 
{vq = O^vs = 1), the Andreev levels are shifted by tt 
to give £j = ±Ay^1 — Vs cos^ ^. The corresponding 
charge current is then 



e VsAsincp 

-- — tanh(— - . 
hy/l-Vs cos2 ((p/2) ^2T^ 



(40) 



The cross-over from the to tt state occurs at vq = vs 
where the Andreev levels are ej = ±A/^/3. These An- 
dreev levels are independent of the phase difference be- 
tween the two superconductors leading to a zero Andreev 
current consistent with equation (37). 

In figures [3] (c) and (d), the current kernel for the left 
side of the interface j^^ (e, cp) is plotted for 'Uq = 1, '^5 = 
and Vq = O^vs = 1, respectively. The current ker- 
nel j^'^ (£,(/?) is integrated over energy £ to give the 
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total current for a specified phase difference, = 
I ^Ja^i^^^)- The kernel j^'^{e^(f) indicates which 
states, i.e. Andreev levels and continuum states, are par- 
ticipating in transporting current through the junction. 
The direction of the current is given by the sign of j^'^. 
As can be seen in panel (c), the lower of the two Andreev 
levels is occupied which is consistent with quasiparticle 
states below the Fermi surface being occupied. The same 
is true for in panel (d), although the current kernel 
has been shifted by tt. 



2. TT junction with a small tilt angle 

For a junction with zero tilt angle, the solution of the 
boundary condition problem reduces to the static spin 
case. The t-matrix equation is simply 

tiit^t') = tiit^t') + [f^^^pJ(^,to (41) 

with f ^ = f = [v^^ og^o v%t,t'). Since the hopping 
elements are time independent, the t-matrix equation can 
be transformed into energy space where the solution can 
be found as 

t°(£) = [[i-f°o5°]-^of°](e). (42) 
The t matrices are then used to calculate the incom- 



ing and outgoing propagators according to equation (16). 
The poles of these propagators subsequently give the An- 
dreev levels, ^j, which turn out to be the same as for a 
static spin {ujl = 0) with arbitrary tilt angle. 

If the classical spin acquires a small tilt angle "0 = 5 
such that sin (5 (5 and precesses around the axis with 
frequency uol^ the spin-dependent hopping element may 
be approximated with vses{t)-a = vs{crz-\-S e~'^^^^^^ ax) 
assuming that vq = 0. This means that a quasiparti- 
cle may either tunnel across the junction with its spin 
and energy conserved, or it may gain or lose energy 
ujl while having its spin reversed as illustrated in fig- 
ure [2] (b). Then, to first order in 5, the t matrices 
are = ^Jta while f = f ^ + ST a- Note that 
STa = e~*^^^jf ^ + e*^^^jf ^ and has no diagonal time- 
independent term. Inserting this into the t-matrix equa- 
tion and disregarding terms of order 6'^ and higher gives 
an equation for Si^ as 

6io.{t,t') = [Staoii+gloii)]{t,t') (43) 
+[tyglo5Q{t,t'). 

Since there are no energy shifts in the last term, the t- 
matrix equation can be written according to equation 
(l24|) as 



1-Af 









At 



5tl\ 













(44) 

where 5Ta^^ = 5T\!^o(l^g^ot^^). The time-independent 
change in the t matrix is consequently zero and the total 



t matrix is given by to, = + e-'^^^St^ + e'^^^Sti. In 
other words, a small tilt angle adds two time-dependent 
components to the t matrix. These two components have 
poles at since quasiparticles gain or lose energy 

ul while being spin flipped. However, these sidebands do 
not show up in the Andreev levels for the charge current 
in this first-order approximation since t^ = When 
the tilt angle becomes large enough for higher-order pro- 
cesses to become important, the sidebands do appear in 
the Andreev level spectrum which is shown explicitly in 
figure [4j Spin-down (spin- up) quasiparticles in the state 
are scattered into the upper (lower) sideband, which 
has an opposite spin direction and the quasiparticles ex- 
change the energy cjl in this spin- flip process with the 
precessing spin. A cross-section of the spin-resolved den- 
sity of states for (/9 = is shown in figure [4]^a). 

In addition to the sidebands, the scattering processes 
also lead to a nonequilibrium population of states. The 
quasiparticles occupying the lower Andreev level £j are 
scattered into states with energies larger than the Fermi 
energy, leading to an occupation of states above the 
Fermi surface. Quasiparticles are also scattered into the 
continuum states below the gap edge. The continuum 
scattering leads to a reduced life-time for the quasipar- 
ticles in the state and can be seen as a broadening 
of the state for energies < — A + cj^A. Despite the 
time-dependent dynamics of the rotating spin and the 
nonequilibrium population of states, the charge current 
is still time independent consistent with reference |22j. 



3. TT junction with arbitrary tilt angle 

As the tilt angle increases, scattering into the side- 
bands increases. The spin-degenerate state £j now splits 
up into a spin-up and a spin-down state and the spin-flip 
scattering still occurs between states separated by energy 
The population of the sharp states, the Andreev 
levels, can be understood on the basis of the spin-flip 
scattering processes. As before, the spin-flip scattering 
connects two states in opposite spin-bands separated by 
energy ujl. This connection results in a population of 
the upper state under the condition that the lower state 
is occupied. In other words, if the spin-down propaga- 
tor is occupied, the spin-up propagator is occupied as 
well. Vice versa, if the spin-down propagator is not occu- 
pied, the spin-up propagator is not occupied either. This 
non-equilibrium occupation is an effect of the precessing 
classical spin which scatters between the two spin states. 
The fine details of this are shown in figure |4] where we 
plot j^'^{s] cp) for some values of the tilt angle 1^. 

This nonequilibrium occupation of states generates a 
current-phase relation quite different from the non-spin- 
flip current-phase relation as shown in figure [5] In the 
laboratory frame, the tilt angle between the precess- 
ing spin and the external magnetic field determines the 
splitting of the Andreev levels. The splitting in turn de- 
termines the current-phase relation and the locations of 
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FIG. 4: (Color online) The charge current-phase relation, 
j^{(p), charge current kernel, j^'^{e,(f), and density of states 
are plotted for a tt junction (vo = 0,vs = I) with tilt angle (a) 
= 0.l7r/2, (b) ^ = 7r/4 and (c) ^ = 7t/2. The temperature 
of the leads is put to 10~^A. The spectral charge current, 
j^'^{s, (p), shows the Andreev levels and their population. At 
some phase differences cpc < ^ < ^^tt — cpc, scattering between 
the Andreev levels and continuum states, which are indicated 
with the black dashed line, cause a broadening of the Andreev 
levels. The charge current, j^{(p) (plotted in units of eA/h), 
is the energy- integrated spectral current and displays abrupt 
jumps at phase differences where Andreev levels become pop- 
ulated/unpopulated. The density of states (DoS) on the right 
sides are plotted for the phase difference cp = and shows the 
splitting of the spin-up and spin-down Andreev levels as well 
as the scattering of the continuum levels into the gap. The 
precessing spin has a frequency of cjl = 0.5A in all plots. 




FIG. 5: (Color online) (a) Current-phase relation for junction 
with {vo = 0,^5 = 1), OJ = 0.5A and i9 = 7r/4. (b) Current- 
phase relations for (vq = 0,^;^ = 1), uul = 0.5A, and varies 
between and 7r/4. (c) cl; = 0.5 A, = 7r/4, t'o = 0, and vs 
varies between 1.0 and 0.3 (d) (vo = 0, = 1), = 7r/4 and 
(jJl varies between 0.00 A and 2. 50 A. As seen sharp step-like 
features in the current-phase relation may be washed out by 
either having a limited transparency t's <C 1 or by having a 
large frequency cjl ^ A. In all four panels the temperature 
is set to zero and the currents are plotted in units of eA/h. 



the abrupt jumps as can be seen in panel |5](b). If the 
classical spin is aligned with the magnetic field, there 
are no spin-flip scattering processes and the only tunnel- 
ing processes taking place are the usual ones where the 
quasiparticles' energies are conserved (as in figure |2] (a)). 
Hence there is no splitting of Andreev levels and there are 
no abrupt jumps in the current-phase relation. If, on the 
other hand, the classical spin precesses in the plane, there 
are only spin-flip scattering processes present. Hence, the 
Andreev levels split up into a single spin-up and a spin- 
down levels. Since each spin band only has one Andreev 
level, there is only one jump in the current-phase relation. 



When the spin-dependent hopping strength, v^, is de- 
creased, the abrupt jumps in the current-phase relation 
disappear as can be seen in panel|5](c). The reason is that 
the Andreev levels are closer to the gap edges for lower 
transparencies. If they are close enough they even merge 
with the continuum states that have been scattered into 
the gap. Similar effects are created when the precession 
frequency cjl is increased. The continuum states are scat- 
tered further into the gap, removing the sharp Andreev 
levels. This modification of the Andreev levels causes the 
jumps in the current-phase relation to be smoothed out 
(panel [5] (d)). 
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c 



(f) 



— (I) 



(I) 

(5) 



:5 



The corresponding t matrix for the right side is found 
by interchanging L and R. The spin-down component 
is given by substituting ^ind reversing the sign of 

the energy shift, —ujl The form of the t matrix 

indicates that a Green's function at a given energy e on 
one side of the interface is connected to Green's functions 
at energies e ^ujl on the other side of the interface. 

As before, the current is carried by sharp Andreev lev- 
els inside the gap as well as continuum states in regions 
:^ujl around the gap edges ±A as shown in figure [4] (c). 
The Andreev levels are given by 



FIG. 6: (Color online) Scattering processes in a Josephson 
junction coupled to a spin precessing with tilt angle ^ = n /2 
and = 0. The tunneling processes are such that all quasi- 
particles traversing the junction reserve their spins while gain- 
ing or losing energy ujl- 



^J,t = -±AVl + /(a;,(^,^,) (50) 



where 



4. TT junction with tilt angle 1^ = 7r/2 

In the limit it/ 2 in which the classical spin pre- 
cesses in the plane, the only spin-dependent tunneling 
present generates spin-fiip scattering processes causing 
quasiparticles gain o r lose energy 0;^ as shown in figure 
[gI (see also reference J^^'^ . Further assuming that the 
spin-independent tunneling is negligible, i.e. = 0, the 
hopping element is 

(45) 



where 



vs sm ? 



and 



Vs sm ? 
















(46) 



which leads to a f q, where only the diagonal term is 
nonzero and has the form 

= [v'og%^^ov^]{t,t') + [v^og%^^ov^]{t,t'). (47) 

Carrying out the convolutions, we find that [v^^ ^ 
v'^]{t^t') is associated with the spin-up quasiparticles 
while [v-^og^ov^][t,t') is related to the spin-down quasi- 
particles. These matrices result in a t-matrix equation 
given as 

1-Af 










1 



Ad 
^0 







1 




(48) 



which immediately results in the only nonzero t-matrix 
term being t^. These two properties ensure that we can 
divide the t-matrix equation into two sets of equations, 
one for each spin band, and Fourier transform the equa- 
tions into energy space. For the left side, we find 



iLt(e) 



1 



(e - wl) -flLt,oo(£) 



(49) 



8vt 



(1 



COS 



2 



\2AJ 



(51) 



A\2 



4vj 



(1-«4)2V (l + «4)2 



(1 + cosipy 



1 



V2Ay 



In the laboratory frame, the number of Andreev levels 
present when i9 = 7r/2 is thus half the number of states 
in junctions with < < 7r/2. The states lie sym- 
metrically around :^ujl/'2 and disperse with the phase 
difference cp until they touch the gap edges and merge 
with the continuum at (pd^L)- The sharp states reap- 
pear again at 27r — (pd^L)- The population of the states 
is indicated in figure El (c) showing the current kernel 
^^'^(e, cp). Since the number of states is reduced, two for 
the incoming propagator and two for the outgoing prop- 
agator, there is only one abrupt jump, instead of two, in 
the current-phase relation for '^^ = 7r/2 shown in figure [s] 

(b). 



5. Tunnel limit 

For a low transparency junction with vq^vs <C 1, we 
may use the first-order approximation ta{t^ t') = Ta{t^ t') 
and obtain an analytic expression for the charge current. 
The charge current between the two spin-singlet super- 
conductors is time-independent and is given by the diag- 
onal hopping element tj^ = T^. In the absence of an ap- 
plied voltage bias, the current is given by the anomalous 
Green's functions and has the form (at zero temperature) 



forwL < 2A, 



eA 
2h 



for UJL > 2A, 



Vs cos2 - ^Vs siii'dKi^) ) sm^, 



(52) 

where K is the complete elliptic integral of the first kind. 
For an alternative derivation, see reference [26 . The first 
term is due to the processes which do not fiip the spin 
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of the tunneling electrons, og^^^ov . The part of the 
current depending on the precession frequency is entirely 
due to the spin-flip processes described by v^og^^ov^ and 



i)-^ o g^ov^ . Equation (|52|) shows that also in the limit 



of low transparency but arbitrary precession frequency, 
the critical current is reduced by spin-flip scattering^ 
and even reversed in junctions dominated by spin- flip 
scattering^^. 



smce 



/ 



27r 



eTT[f3&[ilie,t),g'^ie)]<]=0. (58) 



As a passing note, the charge current is due to the scalar 
component of the t matrix, 

fait) = Ynj ^27riTr[f3K(e,t),5°(e)]<], (59) 



B. Spin currents 



The charge and spin currents in equations (17) and 



(18) are evaluated as the difference between the incom- 



ing and outgoing propagators in each lead a. Defining 
a current matrix as j^{e^t) = g]^^{£^t) — g^^{e^t) and 
using the expressions for the incoming and outgoing prop- 
agators, equation ([l6|, it is easily found that the current 
matrix is given by the commutator 



(53) 



The t-matrices have the spin properties of equation ( |23| ) 
and can, without any loss of generality, be divided into 
a scalar singlet part = ~ ^ spin-vector 

triplet part ii{e,ujL) = tl{e,ujL) ^ e-'^'^Hl{e,ujL) + 
e'''^Hi{e,ujL) such that 



I 



de 
2^' 



-i£{t-t') 



The scalar and spin-vector parts of the t matrix in 
Nambu-spin space are, e.g. for the retarded (R) matrix 
component. 



and 



ASM 



la 



(54) 



(55) 



where 7^'^/(/>^^ and are 2x2 matrices in spin 

space. The spin-triplet t matrix can be parameterized 
similarly to the Green's function in equation (|9|), 



I ^eTv[fs\ii{e,t),gl{e)]<]=0. (60) 



1. Small tilt angle, i9 <^ 7t/2 

Contrary to the charge current case, there are first- 
order contributions to the spin current if the tilt angle is 
assumed to be small, i.e. sini? ^ "d. In section [III A 2[ 



it was found that for a small angle the change in the 
t matrix St compared to the t matrix for zero tilt angle, 
t^, i.e. ia = Pa + ^^a, was given by Sta = e-'^^^^Sil + 
e'^^^^dtj^. Hence, the triplet t matrix is = 6ta giving a 
spin current matrix of 



jl;< =27Ti[5ii,gl]<. 



(61) 



The z component of the spin current is zero since Si^ = 
(as was found in section III A 2 ) . Instead, the spin current 
consists of the two components 



fait) 



-iojLt 



5jl 



5ji 



(62) 



scaling linearly with <C 7r/2 and where 

'^^"^^ = i/ £27riTr[f3<T[<5tlA(e,t),5°(e)]<]. (63) 

The upper two panels in figure [7| show the current kernels 
^^'^(e, (/:?) and j^'^ (e, and the resulting spin-current- 
phase relation for a small tilt angle, ^ = 0.l7r/2. The 
two scattering states carry spin currents in opposite di- 
rections as opposed to the charge current Andreev levels 
which carry a charge in the same direction, see figure [4] 
(a). 



i<Jy(^^ ■ cr) 



(56) 



2. Arbitrary tilt angle 



although the parameterization using the basis (z,t, i), 



e.g. 



la 



e"'"'7l',t^+ + ^'"'^l';^-, has a 



more straight-forward time dependence. 

From the current matrix, the division of the t matrix 
into scalar and triplet components, and the fact that 
is a scalar, it is clear that the spin current is due to the 
triplet components of the t matrix 



de 
87ri 



(57) 



Increasing the tilt angle modifies only the magnitude 
of the spin current. The scattering into the sidebands 
increases as seen in figure [t] (c) and (d) while the z 
component of the spin current remains zero. Conse- 
quently, the spin polarization still rotates with the pre- 
cession frequency ul and can be written as j^it) = 
jt^^-iuLt ji^^iujLt wjiere takes the form 



\ J ^2^zTr[f3(T[tl/^(.,t),^^(.)]<]. (64) 
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FIG. 7: (Color online) The spin current kernel, (a) j^'^ and 
(b) j^'^, for a TT junction ('^o = 0, Vs — 1) with tilt angle ^ — 
Panels (c) and (d) show j^'^ and j^'^, respectively, 
for tilt angle ^ = 7r/4. The precessing spin has a frequency 
of cjl = 0.5 A and the temperature T = 0. The resulting 
spin-current-phase relation is plotted in units of A on top of 
each panel. As seen the phase-dependent population of the 
scattering states makes the spin-current an almost steplike 
function of phase as was the case for the charge-current-phase 
relation shown in figure^ 



Similarly to the charge carrying density of states (fig- 
ure [4|, the scattering into the continuum states increases 
with increasing tilt angle. Figure [S] shows the integrated 
spin currents, i.e. the current-phase relations for the t 
and I components of the spin currents for a transparent 
junction with vq = {) and Vs = 1.0. The abrupt jumps 
in the spin-current-phase relations result from the loss or 
gain of population of spin scattering states as the phase 
difference ip changes shown in figure [71 

Taking a closer look at the triplet commutator and us- 
ing equation (55), the spin-current matrix is (suppressing 
the index a for the leads) 



(65) 



Dividing the retarded t matrix into an anomalous t ma- 
trix, 







lay 



icry{4>A^) 
and "normal" t matrix. 



(66) 



7^^=f^f In, ). (67) 
(similarly for the Keldysh and advanced components). 



FIG. 8: (Color online) Current-phase relations for the spin- 
current components (i^,i^) for vq = and Vs = 1.0. The 
time-independent j^-component of the spin current is zero for 
all cases studied (full black line). The precession frequency is 
(jj = 0.2 A, the tilt angle i9 = 7v/4 and the temperature T = 0. 
The spin currents are plotted in units of A. 



and doing the same for the lead Green's function 

(.n)ia 



fO,R 
J A 



(68) 



.0 _ / < 



the spin-current matrix is given by 



-27ri[<^Ji,/0]< (69) 



where the first term is a contribution from the normal 
Green's functions and the second term is given by the 
anomalous Green's functions. The last two terms do not 
contribute to the current since their nonzero elements are 
off-diagonal in Nambu-spin space. We will separate the 
spin-current into two contributions, one from the normal 
parts of Green's function and t matrix and one from 
the anomalous parts Thus we write 

3n = \j ^27riTr[f3^[7^(£,t),ff^(£)]<] (70) 
j ^27riTr[f3<T[<^^(£,t),/^ (£)]<]. 



and 



3a 



The spin current contains contributions from both as can 
be seen in figure |9] where and j\ have been plotted 
as functions of the phase difference in panels (a) and (b). 
In panel (c), the normal and anomalous contributions to 
the charge current have been plotted for comparison. As 
the figure shows, the charge current is completely given 
by the corresponding anomalous contributions. 
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FIG. 9: (Color online) Spin-current-phase relations for (a) 
the normal spin current, Jn, (b) the anomalous spin current, 
jA, and (c) the normal and anomalous charge currents, 
and for vq = 0, precession frequency u — 0.2 A, tilt angle 
^ = 7r/4 and temperature T = 0. The spin currents are 
plotted in units of A and the charge currents in units of e A / h. 
All currents are scaled by T>s. 



The spin current can be expressed in terms of the di- 
rection of the rotating spin, as 

f{t) = ^Vsl3HCos^{jH)xS{t) (71) 



S 



^/VoVs(3±S^{t), 



where Vq is defined as in section III A 1 The compo- 
nents j3H and are plotted in figure |1Q1 is finite 
only when we have a mixed scattering, i.e. ioi vq > {). 
Moreover, {3^ depends quadratically on the precession 
frequency and has a sin (/?-like dependence on phase and 
is zero for both and tt junctions when in their zero- 
current or ground state. For high hopping-dependent 
strengths, Vs = 1 and = 0, the component has 
a current-phase relation with abrupt jumps which are 
smoothed out for Vs < 1 and > 0- In contrast to 
the component depends linearly on uol for small to 
intermediate precession frequencies, ujl ^ A/2. 



3. Tunnel limit 



In the tunnel limit, the triplet t matrix is given by 
= e-^^^^f ^ + e^^^^f i + f ^. The normal part of the 



spin current is given by the normal parts of the t matrix 
as in equation (70), such that j% = e-'^^^^jj^+e-^'^^^j]^. 
Performing the integral in equation (64) assuming that 
cjl ^ A results in 



^ ^'^1 sini?cosi?|47rzA 



■8iAE 



2A / 



(72) 
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FIG. 10: (Color online) The phase relation of the functions 
Ph and I3± at zero temperature for various realizations of the 
tunneling strengths. The precession frequency is u = 0.2 A 
and the tilt angle is -i^ = n/4. The part of spin current that ro- 
tates with the spin precession ^ /Sh is always non-zero with a 
phase-dependent modulation around a mean amplitude. Note 
the different frequency scalings of 13 ^l) and /3^(^ ^'l)- 

where E{x) is a complete elliptic integral of the second 
kind. A Taylor expansion to second order gives j-^ ~ 

'^1 sin 1? cos 1? . Similar integrations yield = —j^ 
and jz = 0. The normal spin current at zero temperature 
is then 



AS (.^ _ 1 ,.2 ^^L 



cos^(7ff) X S{t). 



S ^16A 

Similarly, the anomalous integrals produce 

= (^—ivqVs sin "dsiiiif^ sin cos cos (p) 



(73) 



8iAK (^) - 47riA 



(74) 



where K{x) is a complete elliptic integral of the first 
kind. In the limit cj^ <C A, the ul dependence reduces 
to 8iAK{^) - 47rzA ^ 7riA{ujL/2A)'^ . Also for the 
anomalous spin current, = and the total anomalous 
spin current is 

m) = -^vl^^^^cosniH)xS{t) (75) 



1 ujj sin u) ^ , , 
To summarize, the total spin current is 



(76) 



fit) = ^VsI3h cos 19 (jH) X S{t) - ^V'DoT)sl3±S^{t) 

(77) 

where (3h = giAp " cos(^]wl, /3_l = Ha^l ^^'^ 
transmission coefficients reduce to Vs ~ 4w| and Vq « 
4^0 for Wo, Vs ^ 1. 
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FIG. 11: (Color online) Temperature dependence of the func- 
tions a, Ph and P± for precession frequency cjl = 0.2At=o 
and tilt angle i9 = 7r/4. The functions are plotted with 
phase difference = tt or = depending on if the junc- 
tion is in a TT state or a state. Note the strong response 
at low temperatures captured in 13 h- At these tempera- 
tures the quasiparticle-caused Gilbert damping is frozen out, 
a{T 0) = 0, and the Andreev- level dynamics alone will 
couple back to the precessing spin. The abrupt jumps at 
T < Tc are when ujl = A(T)/2 and the spin scattering 
connects the square-root singularities in the superconducting 
density of states at s = ^(^)- 

4. Temperature dependence 

The Josephson spin current described above is an effect 
of quasiparticles interfering constructively along closed 
loops when the quasiparticles are subjected to Andreev 
reflection and transmission across the junction as spec- 
ified by the rotating classical spin. The singlet-pairing 
nature of the leads does not support spin currents. This 
means that the spin-dynamics of the Andreev levels will 
be restricted to a small volume (or area) with radius given 
by the coherence length in vicinity of the nanomagnet 
junction. 

However, there is another source of spin current not 
dependent on the Andreev scattering proceses. At tem- 
peratures T > 0, the normal contribution to the spin 
current, given by equation (70), results in an extra 
term 



1 

52 



VsaS{t) X S{t) 



(78) 



which is a spin current carried by quasiparticles and is 
therefore independent of the superconducting phase dif- 
ference if. The parameter a depends in general on the 
precession frequency uol and the temperature T. This 
contribution to the spin current, j^^, may have a finite 
spin-polarized dc component along the z axis. This situ- 
ation is similar to a ferromagnetic quantum dot or layer 
coupled to normally conducting leads where a precession 
of the magnetization leads to the junction behaving as a 



spin pump^^'^-*^'. In the tunnel limit, it can be shown that 
the parameter a = l/27r. 

In figure [11] the temperature dependence of the func- 
tion a(T, ujl) related to quasiparticles is plotted as well as 
the function /3h{T^ujl) which is related to the Andreev- 
level dynamics. At temperatures above the critical tem- 
perature, T > Tc^ the spin current is completely given by 
the quasiparticle spin current j^^ and a = l/27r. As the 
temperature decreases, the normal quasiparticles freeze 
out as the superconducting gap opens and a ^ as 
T/Tc 0. The functions a{T,ujL) and Ph{T^ojl) de- 
pend on the hopping amplitudes vs^ vq as well as on the 
junction state, i.e. the superconducting phase difference. 
For a TT junction, whose tunneling is dominated by vs^ 
a{T^ujL) and Ph{T^ujl) are only weakly dependent on vs 
and vq. The function P±{T^ujl) has a sinusoidal (p de- 
pendence and since it is zero for both tt and junctions 
it is not shown in figure 
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The parameter Ph(T^ujl) 
is zero above the critical temperature and increases as 
T/Tc 0, saturating at ~ f^L<^T>Tc- For junc- 
tions, where vq < vs, the reduction of a(T, cjl) is slower. 
At the same time, Ph{T^ojl) saturates at a lower value, 
^ Yq^lo^t>Tc' The parameter a(T, cjl) does not depend 
on the precession frequency for cjl ^ A/2. 



C. Triplet correlations 



In section IIIB[ it was shown that a rotating classi- 
cal spin inside a phase-biased Josephson junction pro- 
duces a time-dependent spin current. This is somewhat 
surprising since the Josephson junction was assumed to 
consist of two superconducting leads with s-wave symme- 
try. The Andreev processes depicted in figure [2] (b) pro- 
duced by the rotating spin lead to new spin-pairing corre- 
lations which are formed when these scattering processes 
result in positive interference along closed loops. The 
additional pairing correlations to the usual spin-singlet 
ones ~ ^{ip^ip^ — i^ii^^) are the spin-triplet components 

U^t^i^^i^t)^ i^t^t) i^i^i)' The induced triplet 
correlations can be expected to form near the junction 
due to the spin mixing and locally broken spin-rotation 
symmetry provided by S{t)^^^. A similar situation ex- 
ists in SFS junctions with conical ferromagnets; spin cur- 
rents arise due to spin-triplet correlations induced by a 
helical rotation of the magnetization direction in the fer- 
romagnetic layei^^. The spin-triplet correlations induced 
by the rotating spin are localized near the junction and 
are evanescent on length scales on the order of the super- 
conducting coherence length, ^0 = hvf/27TTc. The for- 
mation of triplet Cooper pairs depends on the details of 
the scattering off and the tunneling over the precessing 
spin, such as the precession frequency, cjl, the tilt an- 
gle, and the relative amplitude of hopping strenghts, 
vq^vs- It also depends on the leads through the super- 
conducting phase difference and the temperature, T. 

We now want to quantify the pairing correlations gen- 
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FIG. 12: The components of the d vector, (d^, d^, dz), as functions of temperature. The angle 19 between H and ^S* is ^ and the 
precession frequency is 0.2At=o- In each panel, the spin-independent hopping amplitude Vo is fixed while the spin-transmission 
coefficient T>s is varied. The components are scaled by T>sOJl- In the left panel all junctions are vr-junctions while in the center 
and right panel junctions with Vs = 0.1 and Vs = 0.1,0.2 respectively are 0-junctions. For the 0-junctions the corresponding 
d- vector measures even-in k odd-in e spin-triplet correlations 



erated in vicinity of the precessing spin. Turning first 
to the spin-singlet components we extract the anoma- 
lous Green's functions are ff{±k,e) from the matrix 
[^zn,. ^^ont,.]< ^ 2[gl o (i + [ti,gl]o)]< and write 



de 



[ff{ke)^ff{-ke)] (79) 



with "+(—)" referring to the anomalous scalar matrix 
component of the incoming (outgoing) propagator on 
the left side of the junction and vice versa for the right 
side. V^(/c) is a measure of the (singlet) paring corre- 
lations available to form a singlet order parameter as 

A,(^) = A,7?(^)(7/(^0V^(^0)fc^.n>O' where r]{k) = r]{-k) 
are basis functions of even parity on which the pairing 
interaction may be expanded. The energy Sc is the usual 
cut-off that appears in the BCS gap equation (12). 



The anomalous component of the surface propagator 
also has spin- vector parts, fl^ and /^^^, resulting from 
the spin-scattering processes in the junction. These in- 
duced anomalous components can be quantified in terms 
of a d vector, which is customary to define in relation 
to triplet correlations. In general, a 2 x 2 triplet order 
parameter is given by = d{k) • aiciy and the d vector 
points along the direction of zero spin projection of the 
Cooper pairs^^. For tt junctions, we define a vector do 
which is odd in momentum and even in energy as 



de 



[ff{ke)-ff{-ke)] (80) 



where the direction of the surface normal is n and, for the 
left side of the interface, ff{^k, s) refers to the incoming 



(outgoing) propagator. For junctions, we instead define 
a d vector which is even in momentum and odd in energy, 
dp , as 



d. 



de 



,[/f(^,£) + /f(-^,e)] (81) 



where is the sign of the energy e. This definition is 
based on the anomalous correlations of [^*^'^ ;^out,t]< 



which is related to the matrix in equation (|53|) by 



2[^io + Spin-triplet 
pairing that is even-in k and odd-in s was first considered 
as a candidate pairing state for ^He and is in principle 
not forbidden by symmetr}!^ although not realized for 
superfluid ^He. 

The triplet correlations span the spin space in such a 

way that /f - |(V^tV^i + V^iV^t)^ /f ^ (V^tV^t) and /f - 
Moreover, the instantaneous spin direction of 
the triplet correlations depends on the rotating spin 
leading to a time dependence for the d vector of the form 

^^^^ + d;e^^^^ 



d{t) = dz -\- d^e 



(82) 



The components of the d vector are plotted in figure [T2j 
As expected, the components are finite when the leads are 
superconducting (T/Tc < 1). Setting Vq = ^ leads to the 
magnitudes of the components being equal except for a 
scaling of VqUOl- These properties are modified for finite 
Vo] as the spin-independent tunneling is increased, an 
asymmetry between and d^ emerges and the universal 
scaling disappears. For low temperatures, T/Tc ^ 0.1, 
one can express the d vector in terms of the direction of 
the rotating spin S as 

d(t) = SLS{t) X S{t) + Sh{iH) X S{t) + d,S,. (83) 
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For TT junctions, the term 5z is zero, while 5h = ^ 
for junctions. In the tunnel limit, Vq^Vs <C 1, it 
is possible to find analytical expressions for these S- 
dependent components. For the odd d vector in the 
tunnel limit at low temperatures, 5l^o = 7rVsSm{(p/2) 
and Sh,o = 4:7riVoVsSm{(p/2). The even-in momentum d 
vector, on the other hand, is cut-off dependent and di- 
verges logarithmically with Sc. For the plots in figure [T2j 
we used Sc = 20 A. Furthermore, the relation between d 
vectors on either side of the interface is dR{t) = —dL{t). 



D. Back-action on the precessing spin 

The spin currents on either side of the interface are re- 
lated by j%{t, if) = —jL{t^ ^ difference leading to a 
torque, T{t) = 3\{t) —3%^{t)^ exerted on the rotating spin 
S. The Josephson spin current consists of two spin- vector 
components, jlj^ = jH,L/RilH) x S ^j^^l/rS^. The 
perpendicular spin currents on either side of the interface, 
j±,LS± and j±^RS±^ are equal and therefore cancel. The 
other spin currents, jH,L{"fH) x S and jH^Ri^H) x are 
equal in magnitude but carry spin angular momentum 
in opposite directions, leading to a torque, here called 
Andreev torque, which is given for a single conduction 
channel by 

TAit) = '^'DsPh cos d (jH) X S{t). (84) 

The torque is parallel to the one generated by the exter- 
nal magnetic field H and hence leads to a shift in the 
precession frequency, ujl + ^VsPh cos'd]. As 



was seen in section [ill C , the rotating spin induces local 
spin-triplet correlations near the junction interface and 
the spin-triplet correlations allows the superconducting 
leads to support a spin current even at low temperatures 
when the quasiparticles are frozen out. However, the spin 
current is nothing but transport of spin-angular momen- 
tum and the non-conservation of the spin current results 
in a torque acting on the rotating spin. The shift in pre- 
cession frequency of the rotating spin is therefore a direct 
consequence of the induced spin-triplet correlations. 

The spin current carried by normal quasiparticles 
transport angular momentum from the rotating spin 
into the leads resulting in a damping of the spin 
precessioiJSIHil. This process is the main contribution 
to the Gilbert damping, which has been studied quite 
extensively (see [28] and references therein). As there are 
many possible contributions to the Gilbert damping, it is 
often entered as a phenomenological parametei^^. Here, 
the quasiparticle torque, r^p, is given by the spin current 
in equation ( 78 ) as 



',p{t) = ^aVsS{t) X S{t) 



(85) 



for one conduction channel. Since the torque is perpen- 
dicular to S{t) as well as S{t)^ it leads to an alignment of 
S with the effective magnetic field H. At temperatures 



above the critical temperature, a = (27r)~^, but as the 
temperature decreases, the quasiparticles freeze out and 
the quasiparticle spin current as well as the quasiparticle 
torque vanish as T/Tc 0. This reduction of Gilbert 
damping due to superconducting phase transitions was 
investigated in references |3Q|31j where the Gilbert damp- 
ing in domains with a precessing magnetization was mea- 
sured at temperatures around Tc. 

We propose that the Andreev torque can be mea- 
sured in the following way: Assume that a nanoparticle 
doped with a few magnetic atoms giving it large spin, 
say S ~ 50fi, is placed between two superconducting 
aluminum (Al) leads. The superconducting coherence 
length, ^0 2 /im for Al), limits the area within which 
the time-dependent spin dynamics affects the Andreev 
refiection to ~ tt^q- ^ contact width ~ 100 nm results in 
n ~ 600 conduction channels. The superconducting gap 
A ^ 200 /ieV in bulk Al, but can be made considerably 
smaller in the point contact. The changes in the power 
absorption spectrum of an FMR experiment can now 
be calculated. The width of the resonance peak consti- 
tutes of inhomogeneous broadening, which is due to sam- 
ple imperfections such as anisotropy fields, and homoge- 
nous broadening Ai^hom, which is due to Gilbert damp- 
ing. The homogenous broadening may be expressed aP^ 
Ai^hom = ^Hac where H = \H\ and ac = ^naVg 
is the Gilbert constant!^. For a typical magnetic field, 
H = 60 mT, the macrospin precession is ~ 1 /aeV corre- 
sponding to a frequency of 10 GHz. For a junction with 
Vs ~ 0.1, the freezing out of normal quasiparticles de- 
creases the width of the resonance peak and results in 
a difference in homogeneous broadening on the order of 
Ai/hom(T/T, > 1) - Ai7hom(T/T, ^ 0) - 26 mT. As 
the temperature is lowered, a shift in the resonance peak 
i^o, corresponding to the shift in precession frequency, 
appears due to the Andreev torque. The shift can be re- 
lated to the Gilbert constant as Aujl/ujl = acPH 
Choosing = 7r/4, fSn ^ ^^le low transparency 

limit for temperatures T/Tc < 1/4 which leads to a 
change in the resonance peak of AHq/Hq ~ 2%. This 
ratio could be made higher by increasing the ratio hn/S 
or by increasing the junction transparency. 



E. Magnetization 

The proximity effect is the penetration of supercon- 
ducting order into a normal metaP^ or a ferromagnetP^. 
The reverse case, the so-called inverse proximity effect, 
occurs when ferromagnetic order penetrates into the 
super conduct oi'-^^^. In the present case, where a pre- 
cessing classical spin sits between two superconducting 
leads, there is an induced magnetization in the leads be- 
sides the induced triplet correlations. In the case of a 
static spin, only the spin-singlet correlations and pos- 
sibly the spin-triplet correlations with spin projection 
along the direction of the classical spin contribute to the 
magnetization^^. For a finite precession frequency, also 
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FIG. 13: (Color online) The phase relation of the functions 
/iL ans IJ.H at zero temperature. The precession frequency is 
cj = 0.2 A and the tilt angle is 19 = 7r/4. The function /jh is 
scaled by Vs. 



the time-dependent spin-triplet correlations contribute to 
the magnetization, which can be written on the form 

Mc, = + M^e-^^^^ + Mie'^^K (86) 

The magnetization is given b3^^ 

Ma{t) = ^ I ^Tr[fM9'a<{e,t)+g:'<{e,m, (87) 

where /i^ is the Bohr magneton. From numerical inves- 
tigations it can be found that 



ML{t) 



VsIJLlS X S- 



s 



iJiHinH) X s 



where ii^^ jhh and jnz are functions of temperature, phase 
difference and the hopping amplitudes vs and vq. The 
low-temperature phase relation of the parameters /i l and 



/j^H are plotted in figure [T3| As low temperatures /iz is 
zero. The parameter jhh has a sinusoidal phase relation 
while jj^L shows a cosine behavior. In the tunnel limit, 
/iL = ^(2 - cos (/:)), fiH = ^vovss'mip while fiz = 0. 
The temperature dependence of the parameters /i^ and 
jiiz for a fixed phase difference is shown in figure [T4j The 
phase difference is (/p = tt if the junction is in a tt state 
{vs > Vq) or a state {vs < vq). For these phase dif- 
ferences, the parameter /j^h (not shown) is zero. The 
parameter /i^, on the other hand, is finite for tempera- 
tures 1/2 < T/Tc < 1. The parameter jul ^ (27r)~^ as 
T ^ for TT junctions while /^l ^ (Gtt)"^ for junctions. 
Note that unlike reference [88] , the magnetization is pro- 
portional to the precession frequency and that M = 
for a static spin. 




FIG. 14: (Color online) The temperature dependence of the 
functions /jl and /x^ at phase differences = tt or = 
depending on the state of the junction {vs > vq or vs < vq). 
The precession frequency and the tilt angle are the same as 
in figure [13] 

IV. CONCLUSIONS 

In this work we have presented a detailed analysis 
of the coupling between the Josephson effect and a 
nanomagnet. The magnetization of the nanomagnet is 
brought into precession by an applied external magnetic 
field. The magnetization of the nanomagnet is modeled 
as a classical spin and the precession provides a time- 
dependent as well as spin- active interface which we han- 
dled by formulating a boundary condition within the qua- 
siclassical theory of superconductivity. The boundary 
condition was formulated in terms of a t-matrix equa- 
tion describi ng the scattering between the two supercon- 
ducting lead d^^ l ^^ . We showed that the t-matrix equa- 
tion could be solved in the laboratory frame by preserv- 
ing the explicit time dependence or in a rotating frame 
where the time dependence has been removed by a uni- 
tary transformation. The rotating frame solution is more 
numerically efficient although it introduces a shift of the 
chemical energies of the spin-up and spin-down bands 
of the leads making the laboratory frame solution more 
suitable for investigations of modifications to the super- 
conducting state of the leads due to the time-dependent 
boundary condition. 

As was already described in reference [26 , the Joseph- 
son charge current is time independent despite the mag- 
netization dynamics. Here, we focus on the changes to 
the superconducting state and we show that this sta- 
tionary solution nevertheless exhibits a nonequilibrium 
population of the density of states due to spin scattering 
between the Andreev levels. The nonequilibrium popu- 
lation of the Andreev levels results in a modified current- 
phase relation and causes abrupt jumps as Andreev lev- 
els become populated or unpopulated. The spin current, 
on the other hand, lacks time-independent components. 
Instead, the spin polarization of the current precesses 
around the direction of the external magnetic field. 

We have shown that the spin current is due to an 
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Andreev- level dynamics with induced spin-triplet compo- 
nents. The spin-triplet correlations are due to the break- 
ing of spin-rotation symmetry caused by the precessing 
classical spin. In more detail, the spin-triplet correlations 
are a result of the combination of Andreev retroreflection 
of electron- and hole-like quasiparticles at the supercon- 
ductor surfaces and the spin scattering processes caused 
by the precessing spin in which the quasiparticles undergo 
spin flip while exchanging energy with the classical spin. 
The triplet correlations are evidence of a dynamical in- 
verse proximity effect produced by the time-dependent 
spin-active barrier. 

The spin currents generate a torque on the classical 
spin changing its dynamics. This back-action is due 
to Andreev scattering processes and shifts the preces- 
sion frequency of the classical spin and can in princi- 
ple be tuned using the superconducting phase difference 
across the junction. There is also a source of Gilbert 
damping at finite temperatures due to the presence of 



quasiparticles^. The amount of Gilbert damping can be 
controlled by adjusting the temperature since the quasi- 
particles freeze out as T ^ 0. 

We suggest that the effects of the Andreev torque can 
be measured in an FMR experiment similar to that in 
reference [30]. A measured shift in precession frequency 
would be a probe of the Andreev-level dynamics. Con- 
versely, by controlling the Andreev-level dynamics and 
the inverse proximity effect, one may control the mag- 
netization dynamics of a nanomagnet^ which would be 
useful in spintronics applications. 
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